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Abstract. - The effective mass m* , and the Lande g-factor of the uniform 2-D electron sys- 
tem (2DES) are calculated as a function of the spin polarization (, and the density parameter 
r s , using a non-perturbative analytic approach. Our theory is in good accord with the sus- 
ceptibility data for the simple 2DES, and in excellent agreement with the two- valley Si-2DES 
data of Shashkin et al. While g* is enhanced in GaAs, m* is enhanced in Si. The two-valley 
susceptibility is treated within a coupled-mode (coupled-valley) approach. The coupled-valley 
model is confirmed by comparison with the Quantum Monte Carlo results for a 4-component 
2DES. 



The 2-D electron fluid (2DES) exhibits a wealth of intriguing physics, straddling a rich 
phase diagram [1,2]. The phase diagram contains spin-polarized states at sufficiently large r s , 
say ~ 20 — 27. Here r s = (nn)' 1 ^ 2 is the electron-disk radius [3,4] at the density n, in atomic 
units. It is also equal to the value of the coupling constant T = (potential energy) /(kinetic 
energy). The intermediate regime r s ~ 5 — 20 also hosts many ill- understood phenomena 
including the metal- insulator transition (MIT) [5]. Anomalous values (e.g, see [6]), of g* and 
m* have been found. Some experiments suggest that an enhancement of g* is responsible for 
the strong enhancement of m*g*, while results [7] on Si metal-oxide field effect transistors 
(MOSFETs) suggest that it is m*, and not g* which is enhanced. In this study we show that, 
for ideally thin 2-D layers, g* is enhanced in GaAs-like systems, while m* is enhanced in Si-like 
multi- valley systems. The existence of a coupled- valley state follows naturally from the physics 
of the Si system, and here we present a model leading to excellent quantitative agreement with 
experiment, and with Quantum Monte Carlo (QMC) simulations of a 4-component 2DES [8]. 

Fermi liquid- type theories [9] are valid for r s < 1. Such perturbative methods have been 
applied, invoking impurities [10], or charge and spin-density wave effects [11]. On the other 
hand, QMC calculations of m* involve the excited states of the 2DES and are less reliable 
than for the ground state. QMC results up to r s = 5 have been reported [12]. 
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We showed recently that the 2-DES, 3-DES, and dense hydrogen can be studied using a 
mapping to a classical fluid [13-16]. The accuracy of the map was established by comparison 
with QMC and other independent calculations. Here we use this classical map to evaluate m* 
and g* for the low-density 2-DES. The method is best understood within a density- functional 
picture. 

The density-functional perspective.- The Hohenberg-Kohn-Mermin theorem asserts that 
the Hclmholtz free energy F is a minimum at the true density [17]. If n(r) is the true density, 
it obeys the variational equation 5F[n(r)]/ 6n(r) = 0. If the origin of coordinates is on an 
electron, then if n(r) is the density as seen from this electron, it is a pair- density such that 
n(r) = ng{r). Here g(r) is the electron-pair distribution function (PDF). The variational 
condition gives the Kohn-Sham (KS) equation as usual. Then n(r) is obtained via a sum 
over the KS orbital-densities \ipi\ 2 weighted by the Fermi factors /j. If the electrons formed a 
classical system, the variational equation becomes the Boltzmann form for the density: 

n(r)=ne-M v °°^ +v ^ +v °W. (1) 

V cou {r) is the Coulomb interaction between the electron at the origin and the electron located 
at f. Similarly, V p (r) is the Poisson potential at r, and V c (r) is a correlation potential. 
For a classical system the V xc (r) of standard KS theory is replaced by just a correlation 
potential V c (r). In effect, Eq. ^ evaluates the g(r) of the classical fluid. However, the g(r) 
of a classical fluid is accurately given by the hyper-netted chain (HNC) inclusive of a bridge 
function [18]. Thus, the extended HNC equation is a classical KS equation where V c (r) is 
the sum of HNC+bridge diagrams. The construction of the Bridge diagrams for the 2DES is 
given in refs. [15,20]. 

The classical map has no exchange, and fails as T — > 0. We rectify these lacunae as 
follows. In a system without Coulomb interactions, g(r) should reduce to g (r) which is 
known analytically (at T = 0) or numerically. The first step of the mapping is to introduce 
a potential <fiij( r ) (where i, j are spion labels) such that 4>%{ r ) generates g%{r) when used in 
the HNC equation for ideal electrons [19]. This leads to an exact treatment of exchange. 

Electrons at T = have kinetic energy. Hence the classical map of the quantum fluid at 
T = would be at some "quantum temperature" T q . This is determined by requiring the 
correlation energy e c of the classical fluid at T q be equal to the e c of the quantum fluid at 
T = 0. This may be regarded as a "calibration" of the classical fluid to recover the quantum 
exchange-correlation energy in the r s range of interest. Here we use the e c (r s ) given by QMC 
(Tanatar-Ceperley results for the fully spin-polarized 2DES for r s up tp 30 were used in [15]). 
Once T q , which maps the T — quantum fluid to a classical fluid is known, finite-T fluids are 
calculated from classical fluids at the temperature T c f = (T 2 + T 2 ) 1 / 2 , as justified in ref. [14]. 
We have shown [13, 15] that the classical PDFs are in very close agreement with the quantum 
fluid PDFs obtained via QMC. The success of the method(refs. [13-16]) for 2-D and 3-D 
electrons, hydrogen fluids, and for 4-component 2-D electron fluids as judged by comparison 
with QMC data establishes it to be a well controlled, highly reliable method. The PDFs 
are easily used in a coupling-constant integration for the exchange-correlation free energies 
F xc . Our finitc-T method accurately recovers the low-T logarithmic terms in F c which cancel 
with corresponding terms in F x . This method, based on a classical mapping of the quantum 
calculation to an HNC calculation is called CHNC [1,13-15]. 

Evaluation of m* and g* .- The evaluation of the susceptibility enhancement m*g* uses 
the T = results for the exchange-correlation energy e xc (r s , £). This is expressed in terms of 
e lc (r s , 0) and e xc (r s , 1), and a polarization factor P(r s ,() given in Eq. (6) of Ref. [15]. Using 
Hartree units, the ratio of the static spin susceptibility to the ideal (Pauli) spin susceptibility 
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Fig. 1 - The compressibility ratio K° /K calculated from the 4-component (2- valley) QMC results of 
Conti et al., and from the 2-component (1-valley) data using the k — ► limit of iTo/n cm where II cm 
is the coupled-mode polarization function. 

is: 

Xp/Xs = (mV)- 1 = 1 + r 2 s d 2 e xc /8( 2 . (2) 

The effective mass m* at temperature T is the ratio of the specific heats, C V (T)/C®(T) of the 
interacting and non- interacting 2DES. 

m*(T)= C v (r s ,(,T)/C (r s ,(,T). (3) 

The specific heats are obtained as the second-T derivatives of the interacting and ideal 
Hclmholtz F(r s ,T). Here T is the physical temperature and not T c f. The latter is used 
only in the classical map to obtain the PDFs. F x (r s ,T) has a logarithmic term of the form 
T 2 log(T) which is cancelled by a similar term in F c (r s ,T). That is, 

F x = A x + BJ 2 log(i) - C x t 2 , t = T I Ep (4) 
F c = A c + B c t 2 log(t)-C x t 2 , B x = -B c . (5) 

This cancellation holds to 85-95% in our numerical CHNC results, for the range r s = 5 — 30, 
< t < 0.25. Thus, at r s = 15 and 25, (B X ,B C ) arc (-0.0258, 0.0228), and (-0.0155, 0.0142). 
If Hubbard-type finite- T RPA were used in the self-energy, the cancellation is quite poor, 
even at low-r s . These logarithmic terms and the m* have also been studied by Geldart et 
al., using CHNC [21]. Multi-valley systems- Shashkin et al. [7], also [23], have studied clean 
low-density 2-valley 2DES in Si-MOSFETs. The two valleys are assumed degenerate [24]. It 
is found [7] that the m* is strongly enhanced, while g* shows little change. The enhanced m* 
is independent of (. These results, "contrary to normal expectations", are reproduced by our 
coupled-mode theory of two valleys. 
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Two equivalent valleys and two spins imply 10 different PDFs, g™, where u,v are valley 
indices. Such a calculation for each r s ,£, T and many values of the coupling constant is 
laborious. A simpler procedure using just three PDFs is possible. Even if £ ^ 0, each valley 
has a density n/2. Thus the 2-valley system may be made up from the known properties of the 
one-valley (two-spin) 2DESs coupled together by their Coulomb interaction. The individual 

1- valley correlation free energies F", are known from QMC and CHNC results. The inter- 
valley term for a system with a total density n, and valley densities n/2 is not known. Here 
we present a simple approximation validated by calculating the 2-valley compressibility in 
the same way and comparing with the QMC data of Conti et al [8]. There is no exchange 
interaction between up-spin and down-spin electrons in the one-valley system, and the spin 
densities are n/2 at £ = 0. Hence, since F"(n/2), F"(n/2), and F c (n,( = 0, [512]) for the 
one-valley system are known, we build up the 2-valley system within the assumption that 
F c (n,C — 0) can be used for the inter-valley contribution to the F c of the 4-component (i.e., 

2- valley) system. In a full 4-component CHNC calculation, the inter-valley interaction is 
switched on via a coupling constant integration. This effect can be recovered within linear 
response by developing the coupled-mode 2-valley response functions. An analogous coupled- 
mode problem arises in electron-hole systems (see Vashishta et al [22]). 

The total (spin or charge) density-fluctuation spectrum of the electrons in a given (single) 
valley v is described by the response functions \ v = x^/D v , where Xv is the 2-D Lindhard 
function weighted appropriately with the square of the Bohr magneton [Ag or unity, and D v is 
a corresponding denominator for each case. The Pauli susceptibility xp is the long-wavelength 
limit H%Xv(k = 0)- Let us consider a denominator of a response function (which may be the 
charge response x, the proper polarization function n, or the spin susceptibility Xs, depending 
on how the local-field factor G v is specified). The denominator D v = 1 — v cou (l — G v )Xv and 
defines G v , the local- field factor (LFF, see [26]). We are only concerned with the static k — » 
limit. Then G v for n are related to K°/K, while the G v for Xs is given by Xp/Xs, a s in Eq.|2| 
and depends on the correlation free energy F c of the one- valley 2DES. When two such 2DESs, 
described by Xv and Xu interact via the inter-valley term, coupled modes are formed. These 
modes are described by the zeros of a new denominator of the response function of the total 
2-valley system. This coupled-mode form is [22]: 

Xcm = [x°u + X v +v 2 couX°uX°^G uv )]/D cm (6) 
^G uv = G u + G v G uv — G vu (7) 
D cm = D U D V - v 2 cou x° u x° v (l ~ G uv ){\ - G vu ) (8) 

Here G uv is an LFF arising from the inter-valley term F uv already discussed, and modeled 
by F\2(nX = 0) at k = 0. Hence we express the susceptibility enhancement Xs/Xp as 
Xcm/xP-, an d this is evaluated from the G U ,G V and G uv . Equation determines G u = G v , 
where the correlation part involves the second derivative (r 2 sv d 2 F% / dC/ 2 ) . Similarly the cross 
term G uv involves r 2 d 2 Fi 2 (n,( — 0)/d( 2 . The 4-component QMC results of Ref. [8] for 
F c (r s ,( = 0, T = 0) enable us to calculate the compressibility ratio K°/K of the 2-valley 
system directly. The coupled-mode theory, applied to the proper polarization function n 
gives another evaluation K°/K. The agreement between the two methods is shown in Fig.^ 
A similar comparison for Xp/Xs is not possible as the QMC results are available only at 
C = 0. However, the agreement between the two estimates of K° /K validates our coupled- 
mode evaluation of 2-valley properties from the 1-valley energies. Thus the 2-valley results 
are constructed from the 1-valley CHNC energies (which argee closely with QMC data) which 
include the usual bridge contributions [15]. 

Results- In Fig. we show Xs/xp — rn*9* for a single- valley system, as a function of the 
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density n, and as a function of r s (see [3]) at T — for £ = 0. Our results, the experimental 
data of Zhu et al. [6] , and QMC data, extracted from Fig. 2 of ref . [2] are displayed. The high 
density regime [27] is in agreement with standard theories and is not displayed. 

For the Zhu et al. data we use their fitted form m*g* = (2.73 + 3.9 n()n~ 0A where the 
density n is in units of 10 10 cm" 2 . The strong agreement between CHNC and the Zhu data 
is perhaps fortuitous since the results are quite sensitivity to the d 2 /dC, 2 calculation to the 
energy differences AE = E c (( = 1) — E c (( = 0) and the form of the polarization factor 
P(r s ,C). The CHNC is calibrated to the Tenatar-Ceperley QMC which differs somewhat 
from the Attaccalite data. We have also plotted two CHNC curves where (see Eq|2J) the term 
r 2 d 2 e xc /d( 2 has been modified by ± 2%. Clearly, errors in converting to r s , modification of 
exchange-correlation gradients by well- width effects and the presence of impurities etc., can 
produce such a change. The bottom panel (Fig. [2Jl shows the comparison against r s . 

Zhu et al. report a £ dependence, but now they consider that the finiteness of the 2-D layer 
and orbital effects cannot be ignored in analysing field-dependent data [28]. As pointed out 
via the ±2% plots in Fig0 the sensitivity of xp/Xs to small errors in the xc-energy gradient 
is also important. (Discussion of these and other data for m* and g*, of Zhu's thesis [6] will 
await their publication). 

In our results, Xs/xp is l ess sensitive to £ at high density, and very sensitive to £ at low 
density, approaching the para-^jerro transition. In fact, the second derivative in Eq[5] diverges 
at C= 1. 

A very different experimental picture is found in Si-2DESs [7]. The CHNC results for the 
coupled 2-valley 2DES are shown in Fig. EI The top panel compares the m*g* = Xcm/xp 
obtained from experiment and the coupled-mode analysis (the LFFS used are for the spin-spin 
response). The inset shows the shift of the simple uncoupled- valley curve to higher densities 
when the valley coupling is introduced. The conversion between density and r s is discussed 
in ref. [3]. 

The middle panel (fig. GJ) shows the m* calculated from the finite-T analysis, with the 
sharp rise occurring at r s ~ 5.4, i.e., density n cv = lxl0 11 /cm 2 . The inset shows the lack of 
£ dependence in m* for three densities. This is because the physics is dominated by singlet 
interactions, as in the ambi-spin phase reported earlier [I]. The lower panel of Fig. [3] shows 
the flat g* of the coupled-valley fluid, while the inset shows the usual increase of g* in the 
uncoupled system as the density is reduced. 

Conclusion- We have presented results for the effective mass m* , and the Lande g* factor 
of 2-D electron fluids, using an analytic method. Our results suggest that exchange effects 
dominate as r s increases in 1-valley 2D system, enhancing g* when the one- valley spin-response 
diverges. Correlation effects outweigh exchange in 2-valley systems where m* is strongly 
enhanced and only weakly dependent on £. The tendency to form singlets already noted in 
the single 2DES [25] becomes stronger in the 2-valley 2DES where a coupled-valley state is 
formed. Our theoretical results depend only on the e xc (r s ,£,T) used in calculating £ and T 
derivatives, and invoke no fit parameters specific to this problem. However, the inter-valley 
energy F uv was approximated via the -Fi2(£ = 0) of the one- valley inter-spin energy. As 
already noted, x/xa is quite sensitive to the evaluation of d 2 /dC, 2 . However, the agreement 
of the present model with experiment may prove useful in understanding the experimental 
results. 

* * * 

We thank Peter Coleridge, Wally Geldart, Frangois Perrot, Sasha Shashkin, Horst Stormer 
and Jun Zhu for their comments and correspondence. 
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Fig. 2 - The spin-susceptibility enhancement Xs/xp = m *9* m the 2DES. Top panel: comparison of 
experiment [6], QMC [2]. and CHNC. Curves marked ±2% are CHNC predictions if the exchange- 
correlation contribution r 2 d 2 e xc /dQ 2 is modified by ±2%. Bottom panel: CHNC results for m* g* for 
3 spin-polarization (, and the experimental £=0 data, plotted against r s . 
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